ON PERIODIC SOLUTIONS FOR 
A REDUCTION OF BENNEY CHAIN 

MISHA BIALY 

Abstract. We study periodic solutions for a quasi-linear system, which 
00 ' is the so called dispersionless Lax reduction of the Benney moments 

^^ . chain. This question naturally arises in search of integrable Hamilton- 

ian systems of the form H — p 2 /2 + u(q,t) Our main result classifies 

completely periodic solutions for 3 by 3 system. We prove that the only 

5— I ■ periodic solutions have the form of traveling waves, so in particular, the 

*— p' potential u is a function of a linear combination of t and q. This result 

"^^ implies that the there are no nontrivial cases of existence of the fourth 

_i_ ■ power integral of motion for H : if it exists, then it is equal necessarily 

to the square of the quadratic one. Our method uses two new general 

observations. The first is the genuine non-linearity of the maximal and 

minimal eigenvalues for the system. The second observation uses the 

\^ | compatibility conditions of Gibonns-Tsarev in order to give certain ex- 

Xy^ • actness for the system in Riemann invariants. This exactness opens a 

_J ' possibility to apply the Lax analysis of blow up of smooth solutions, 

^3 which usually does not work for systems of higher order. 



1. Introduction 

r ■ Let H = p 2 /2 + u(q, t) be a Hamiltonian of a 1, 5-degrees of freedom sys- 

tem with the potential which is assume throughout this paper to be periodic 
function in both variables. There is a conjecture attributed to G.Birkhoff 
saying that the only integrable plane convex billiards are ellipses. The direct 
analog of this conjecture for the Hamiltonian system with 1,5 degrees of free- 
dom would be the claim that the only integrable Hamiltonian functions of 

(^ , the form H = p 2 /2 + u(q,t) are those having the potential functions u which 

are periodic functions of the form : u = u(mq + nt). There are several at- 
tempts to approach this problem. Let me refer to the works [5].[2].[5].[13j.[I] 
for various approaches to this circle of questions which remain outside the 
discussion of this paper . 

In the present note we restrict this question to the search of the additional 
integrals which are polynomial with respect to the momenta variable p with 
the coefficients which are periodic functions of q and t . More precisely we 
want to find all those potential functions u(q, t) for which there exists an 
additional function F(p, q, t) invariant under the Hamiltonian flow (such an 
F is called the first integral of motion) which is a polynomial in the variable 
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p of a given degree, say n + 1, having all the coefficients periodic in q and t. 
Write 

F(p, q, t) = U- lP n+1 + u p n + u^- 1 + • • • + u n , 
and substitute to the equation of conservation of F. 
(1) F t +pF q -u q F p = 

Equating to zero the coefficients of various powers of p, one easily obtains 
the following information. The coefficient u_i must be a constant, which 
will be normalized to be ^px- Also uo must be a constant, which we shall 
assume to be zero (this can be achieved by a linear change of coordinates 
on the configuration space T 2 ). Moreover the coefficient u\ satisfies (u\)q = 
(u)q. Therefore, u\ and u will be assumed to be equal (the addition of any 
function of t to the potential u does not change the Hamiltonian equations). 
Moreover, the column of the of the rest of the coefficients U = (ui, ■■,u n ) t 
satisfy the following quasi- linear system of equations. 



(2) U t + A(U)U q = 0, A(U) 



( -1 ••• 0\ 

(n-l)ui -1 ••• 



2u n _ 2 ••• -1 
V u n _i • • • / 



Notice that the derivative F p of F coincides with the characteristic polyno- 
mial of -4(17). 

In fact the system ([2]) is very well known among integrable systems com- 
munity. This is the so called dispersionless Lax reduction of the moments 
Benney chain ( see for example [6], [7], [12], [14| and references there in). 
There are many beautiful properties of this reduction. For example, this is a 
Hamiltonian system of Hydrodynamics type (in the sense of Dubrovin and 
Novikov see [9]. Moreover, it has infinitely many additional conservation 
laws \j. The most important property for this paper is that the system ^ 
is diogonalizable, i.e. can be written in the form Riemann invariants (see 
([3]) below). Almost nothing is known, however, about the global existence 
of smooth solutions for this system. In the theory of quasi-linear hyper- 
bolic PDEs it is a well known problem to prove the occurrence of blow up 
of smooth solutions. It was performed first by Lax [11] for 2 by 2 systems 
satisfying the so called genuine nonlinearity condition. His method relies 
heavily on the possibility to write the 2 by 2 system in the diagonal form, 
and also on the genuine non-linearity condition. Lax analysis was performed 
in [3] for 2 by 2 system of the form ([2]) where it was proved that the only 
periodic solutions for that case are constants. For systems of higher size the 
original method by Lax does not apply in general. We will show bellow two 
new observations concerning our quasi-linear system. The first is, that for 
the hyperbolic case, i.e. the case when all eigenvalues of the matrix A(U) 
are distinct and real, it follows that the minimal and maximal eigenvalues of 
A(U) are in fact genuinely non-linear in the sense of Lax (see Corollary 12.21 



I was told by M.Pavlov that this is a well known fact among the specialists in integrable 
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bellow) . The second key fact is that the so called Gibbons-Tsarev compat- 
ibility system (the equation (0} of Lemma 12.11 and Corollary I2.3|) provides 
certain " exactness" of the system (see the equation (|6|) and therefore enables 
one to perform the Lax original analysis for higher values of n. 

Our main application of this approach in this paper is the following clas- 
sification of smooth periodic solutions for the system. 

Theorem 1.1. Let n = 3. Then the only periodic solution of the quasi- 
linear system |||) are the traveling waves solutions, where u\,U2,u^ do not 
depend on t. 

Corollary 1.2. Let F = l/4p 4 + u\p 2 + U2P + U3 be a polynomial of degree 
four with periodic coefficients which satisfies the equation F7J). Then there 
necessarily exists a quadratic integral, i. e. the energy H, and F is a function 
ofH. 

In order to prove this theorem we shall divide between different regions: 
strictly Hyperbolic region f^ where all three eigenvalues of A(U) are real 
and distinct, Elliptic region f2 e , where two eigenvalues are complex conjugate 
and the third one is real, and finally the region of degeneracy Oo, where at 
least two of the eigenvalues collide. It turns out that in all these regions the 
behavior of solutions can be understood completely. The proof of Theorem 
II. H is obtained just by patching together the information of these three cases. 
Organization of the paper is as follows. In the next Section we shall explain 
the two basic observations mentioned above concerning hyperbolic situation. 
In Sections O [U [5] we study the regions Qh,Q e ,Qo respectively for the case 
of n = 3. In Section [6] we prove formulas for derivatives of the eigenvalues 
and verify, for the sake of completeness, the Gibbons-Tsarev compatibility 
conditions. 
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2. Main Observations 

Let me denote by Ai,..,A n the roots of the polynomial F p . And let 
rj = -F(Aj) be the corresponding critical values. The starting point for 
us is a beautiful classical theorem by MacLane, stating that the mapping 
{u\...u n ) >—¥ (n, .., r„) is in fact a global diffeomorphism between the domain 
of strict Hyperbolicity (that is the domain of all U = (ui, . . . , u n )' where all 
the roots of the polynomial F p are real and distinct) with the domain of all 
possible critical values in lR n that is of all those (n, ... ,r n ) such that the 
differences (r^ — r^ + i) have the sign (— l) fc + n for all k = 1, . . . ,re — 1 (we 
refer to [10] for E.B. Vinberg's proof of this theorem and further results and 
discussions). According to this theorem (ri,..,r n ) can be taken as regular 
global coordinates in this domain. 
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The importance of these coordinates for our system follow from the fol- 
lowing computation. Substitute p = \ into the equation ([I]). One gets 
the following diagonal system on the variables Ti (they are called Riemann 
invariants) 

(3) (ri)t + \i(ri,...,r n )(ri) q = 0, i = l,..,n 

The derivatives of the roots A, with respect to the critical values r, satisfy 
the following relations. 

Lemma 2.1. The following formulas hold true 

(a) d r .\i = ~ Fpp{Xi ) Efe=i,fc^ aFaT' 

(b) d r .\j = ~ Fpp \ Xi) Xj -x^ i + 0, 

(c) Fp P (Xi)u r . = I, 

(A) 1, — 2Ur i Ur k 

(UJ U nrk — ( Aj _ Afc )2- 

Remark 1. The condition (JdJ) is in fact the so called Gibbons Tsarev com- 
patibility condition (see [7]) (note there is a missprint in their formula-the 
factor 2 is missing). I didn't find however the formula (jaj) in any paper on 
the subject. We suggest the proof of all of them in Section [6] in a very short 
way. 

Corollary 2.2. In the strictly hyperbolic region Ai < A2 < • • • < X n the 

maximal and minimal eigenvalues are genuinely nonlinear in the sense of 
Lax: 

a ri Ai/0,9 rn A„/0 

Another important consequence of the Gibbons-Tsarev conditions is the 
following 

Corollary 2.3. In the Hyperbolic region introduce the functions 

G i = -2 lo gK*l = 2 log 1-^(^)1 
Then it follows from the lemma that 

G = Urirj = — = j 

rj % 2u n (A, - Aj) 2 Aj - Aj 

In order to perform the blow up analysis for our system we shall differen- 
tiate the quantities Wi,Wi = (rj) g along the integral curves of the family Aj. 
They are, by definition, the integral curves of the equation q + Xi(q,t) = 
on T 2 . Let Vi = (1, Xi(q,t)) be the i-th vector field and let L v . = dt + Xid q 
denotes the Lie derivative along the field V{. Differentiating with respect to 
q the i-th equation of ([3]) one gets the following 

(4) L Vi (wi) +w 2 i {X i ) n +w i ^2(X i )r j (r j ) q = 

Notice that by definition: 

L v Tj = ( rj ) t + Xi{rj) q 
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Subtract from this expression the j-th equation of ([3]) 

= (rj)t + \j(rj) q 
one verifies that 

(Tj)q- : 



Aj — Xj 
Substitution of this expression into ((2|) leads to: 

(5) ^.W+^lAik + Wi^lA,)^ = 

Therefore it follows from the last (Corollary 12. 3D that §5$) can be rewritten 
in the following way 

L Vi (Wi) + wf(Xi) n + Wi yX G i)r-i L Vi r 3 = °- 

Therefore, taking into account the i-th equation of ([3j) we get: 

(6) L Vi (wi) + wj(\i) r . + WiL Vi Gi = 

Multiplying by exp G% this equation one rids of the linear term as follows: 

(7) L Vi ((expGi)( Wi )) + (exp(-G,)(A i )r l )(exp(2G i )u; 4 2 ) = 
Using the explicit expression for Gi of Corollary 12,31 and denoting 

Zi = |i ? pp (Aj)| 1/2 Wi = |F pp (Aj)| 1/2 (ri) q , Ki = \F pp (\i)\~ 1/2 (Ai) n 

we get the following equation, which is crucial for the analysis of the blow 
up of the solution. 

(8) L v .Zi + Kizf = 0, i = l,...,n 

As an immediate consequence of this equation we state the following 

Theorem 2.4. Let U = (ui,...,u n )* be a periodic solution of the quasi- 
linear system $&j) corresponding to the strictly Hyperbolic regime, i.e. all 
eigenvalues are real and distinct: Ai < A2 < • • • < A n . Then the Riemann 
invariants n and r n corresponding to the minimal and maximal eigenvalues 
are constants. 

Proof. This fact follows immediately from the equation (jHJ). Indeed, in the 
region of strict hyperbolicity the F pp (\i) does not vanish and so, by the 
genuine non-linearity of Ai and A n the functions K\,K n are bounded away 
from zero. Then it follows from the explicit formula for the solutions of (jH)) 
that the only solution which does not explode in a finite time is zi,z n = 0. 
Thus (7*1 )g = (r n )q = and the equations Q imply that n and r n must be 
constants. This yields the result. □ 

A refinement of this argument is the content of the next section on the 
Hyperbolic region, for 3 by 3 system. In what follows we shall assume that 
n = 3, i.e. the system is 3 by 3. Let me denote by fl^ be the region of strict 
Hyperbolicity and f2 e be the region, where there are two complex conjugate 
eigenvalues for A(U). The complement, T 2 — ($7^ U £l e ) is the set of those 
points (q, t) where the matrix A(U) has at least two equal eigenvalues . We 
shall denote this set Qq. And finally, Oqo will denote the set of maximal 
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degeneration, i.e. where all three eigenvalues are equal and thus equal to 
zero (since the sum of all the three eigenvalues vanishes) . 

3. Hyperbolic region f^. 

Before stating the main result of this section, let me rewrite the formulas 
of the Lemma EI for the case of 3 by 3 system. 

In this case since Ai + A2 + A3 = we get the following simplifications 

(Aik=_ ^F' Kj= "^i^F' F -^=n^-^) 

For the case n = 3 we have the following refinement of the theorem of the 
previous section. 

Theorem 3.1. Let U = (u\, 1*2,113) be a periodic solution of the system 
|Ip ; and let Qh Q T 2 be the domain of strict Hyperbolicity. Then the Rie- 
mann invariants ri,rs are constants in every connected component of the 
Hyperbolic domain Q^. 

Proof. We give the proof for n, the other case is analogous. The first step 
of the proof is the fact that the derivatives {ri) q are bounded functions on 
the domain 0,^. Indeed, by definition, 

n = F(Xi, q, t) => (n) q = F g (Xi, q, t) 

and thus by the periodicity of the coefficients of the polynomial F, all roots 
Xi of the derivative F p are bounded and so are the derivatives (rj) g . Consider 
the integral curves of the Ai-family in the domain fi/j. Suppose that such 
a curve approaches the boundary of Q,h, then F pp — > while (r*)g stays 
bounded. Therefore, it follows that 

z i = {n) q \F pp {\ i )\ 1 l 2 ^0. 

I claim that then z\ equals zero identically in Q^. If, for example, z\ is 
positive at a point (qo,to) then by the equation (jHJ) and the fact that K\ 
is positive, we have that z\ is a decreasing function of time and thus in 
the backward time along the integral curve cannot approach the boundary, 
because on the boundary z\ vanishes. On the other hand, if the integral 
curve stays inside O/j forever in the backward time, then the function K\ 
stays bounded away from zero, and therefore the blow up of the solution 
must occur in a finite (backward) time. Thus z\ can not be positive. The 
opposite case, when Zi(qo,to) can not be negative is completely analogous. 
This argument in a more precise form looks as follows. Denote by M\, Mi 
positive constants such that 



\{r 1 ) q \<M 1 , ^\F pp {\i)\ < M 2 . 



Then for any backward time along the integral curve the monotonicity of z\ 
implies that 



M 1 J\F~(X)\ > {n) qy f\F~(M\ = Zl > zi(go,to) 

and therefore 

— — — £1(^*0) 
F pp( X i)\ > Mi ■ 
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In addition, Ai can not be too close to zero. Indeed, if < — Ai < a, then 
by the zero sum condition also | A2 1 < a and < A3 < 2a. Then one would 
get 

^\F pp (Xi)\ = ((A 2 - Ai)(A 3 - Ai)) 1 / 2 < V6a. 
So together with the previous estimate this implies that |Ai| > ^,° ■ 



Then 



K _ 3 I Ai| 3^i (gp, t ) 

1 ((A 2 - Ai)(A 3 - Ai))V2 - ^q Mi MI 



So K\ is bounded away from zero, and again by the explicit formula for 
the solution of ([8]) it explodes in a finite backward time. This proves the 
claim that z\ vanishes identically in 0,^. Thus z\, {r\) q = 0, and so by the 
equations ([3]) n must be constants. D 

The next theorem describes completely the solutions of the system in the 
Hyperbolic region Q^ . 

Theorem 3.2. Either the solution U = (ui,U2,us) t is a constant solution 
for (0j on T 2 , or the region £7^, is a union of strips on the torus parallel to 
the t -axes and the following relations hold 

u\ = u\(x),U2 = 0, M3 = u\ + const, 

so that the polynomial F equals (up to a constant) in Vth to the square of 
the Hamiltonian 

(2 \ 2 

— + ui(q)\ + const. 

Proof. Let me note, that since r% are successive critical values of the poly- 
nomial F p , then everywhere in Q^ holds r 2 > r\,r<i > r%. On the boundary 
dSlh two of the eigenvalues collide, say A 2 collides with Ai (or with A3 ,or 
maybe both), and hence r 2 = r\ ( or r 2 = r 3 ). By the previous theorem 
ri,r3 are constants in f^, and in addition r 2 has constant values along the 
integral curves of the A 2 -family. Then it follows, that non of these curves 
can approach the boundary, because otherwise this would give r\ = r 2 in 
the inner points. Moreover, the function r 2 satisfies the equation 

(r 2 )t + A 2 (r 2 )(r 2 ) g = 0, 

where A 2 (r 2 ) depends only on r 2 , since n, r$ are constants in Q^. Therefore, 
the characteristics of this equation are the straight lines in the (q, £)-plane, 
and thus there are two possible cases. The first case is that there exist 
two intersecting straight lines of the family. Then r 2 has to be constant 
everywhere ( since r 2 has constant values along characteristics, and any 
straight line intersects at least one of the two intersecting characteristics). 
So in this case all ri,r 2 ,r3 are constants everywhere in fi^, and hence also 
the coefficients of the polynomial u\, u 2 , U3 are constants. In such a case Cl^ 
must the whole T 2 . 

In the second case all the characteristic straight lines are parallel with 
the same A 2 = [i = const , and thus the domain O/j in this case is a union 
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of parallel strips with the slope \x. Next we claim that if the solution U = 
(til, u 2 , Us) f is not constant in £\, then fi = 0. To see this, let me recall that 

F p (fi) = fi 3 + 2uifi + u 2 = 0, 

and therefore 

(9) u 2 = —2[iu\ — /j, 

Substituting @ into the first equation of the system ([2]) we have 

(Ul)t = -(U 2 )q = 2/i(ul)g 

In addition we have that along any characteristic straight line of the family 
A2 the values of ri,r2,r3 are constants, and then also Ai, A2, \3,ui,U2,Ua, 
because ri, r 2 , r% are genuine coordinates. In order to prove the claim let us 
assume that on the contrary \x 7^ 0. Then u\ has to be globally constant in 
fi/i because it satisfies the following two equations 

(10) («i) t - 2/x(«i), = 0, (m) t + \x{ux\ = 

But then, by ([!]), also «2 is constant in 0,^. Therefore Ai,A2,A3 are all 
constants, since the polynomial F p has constant coefficients. Also since 
113 = n — -F(Ai) then also U3 must be a constant. So the solution is in fact 
a constant solution contradicting the assumption of the claim. 

Thus, we get that /i = 0. In this case all characteristics of the second 
family are parallel to the t-axes and the region $7^ is the union of strips 
parallel to the i-axes. Moreover by (jlOD u±,U3 have to be the functions on 
q only and by ([S]) u 2 = 0. 

Therefore F p = p(p 2 + 2u\), and so 

1 fp 2 \ 2 

F = -p 4 + u\p 2 + u 3 = [ — + ui) +u 3 -uf, 



and since ri,r 3 are constants, then u% — u 2 = constant, and we are done. 
Notice that on the boudary of f^, u\ vanishes and so d£lh C £1$q. D 

4. Elliptic region VL e 

Consider now the region £l e where the polynomial F p has two complex 
conjugate roots, say Ai,2 = a ± i[3 with j3 > in Q e , and A3 E R. In this 
case ri t 2 are also complex conjugate, say n t2 = v±iw and r^ is real. Notice, 
that for the points of the boundary of £l e we have Ai = A2 = — A3/2 are real, 
ri = r2 are real also, and so j3 = 0, to = 0. For the region Vt e we have the 
same description of solutions as in the strictly Hyperbolic domain, but for 
completely different reasons. 

Theorem 4.1. Either U = (ux,u 2 ,Us) t is a constant solution for the system 
HP on the whole T 2 , or the region VL e is a union of strips parallel to the 
t-axes on the torus and and the following relations hold 

u\ = ui(q),u 2 = 0,v,3 = n 1 + const 

So that the polynomial F equals (up to a constant) in Q e to the square of 

the Hamiltonian 

2 \ 2 

p \ 

— + ui(q)\ + const. 
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Proof. The Riemann invariants 7-1,2 = v±iw satisfy the equations ([3]), which 
are equivalent to the following elliptic system on their real and imaginary 
parts. 

v t + av q - /3w q = 
wt + (3v q + aw q = 0. 

This is an elliptic system, and therefore, since w vanishes on the boundary, 
then by the strong maximum principle the function w must vanish identically 
in the whole Q e . Substituting back to the elliptic system we get that v is a 
constant everywhere O e . Therefore, Ai 2 are roots of the polynomials F — v 
and of F p . Then, 

F=\(p- Ax) 2 (p - A 2 ) 2 + v = i((p - a) 2 + (3 2 )) 2 + v. 

Moreover, a vanishes because F does not contain the cubic terms. Therefore 

F = (p 2 /2 + (3 2 /2) 2 + v, A li2 = ±i/3,A 3 = 

This situation can be completely analyzed, because in this case the quadratic 
function F = p 2 /2 + (3 2 /2 is the conserved quantity. But then, by the 
equation ([1]) for F we have: 

((3 2 )t = 
((3 2 /2- Ul ) q = 0. 

Notice, that on the boundary of Vt e we have /3 = u\ = «2 = 0- By the first 
equation (3 vanishes on any line q = const whenever it crosses the boundary, 
and it is a constant on any line which lies entirely inside f2 e . This yields 
immediately that f2 e is in fact union of strips parallel to the t axes and 
(3 = (3{q) and also j3 2 /2 — u\ = because on the boundary both and u\ 
vanish. So we proved 

F = (p 2 /2 + ui{q)) 2 + const 

and this completes the proof □ 

5. Degenerate regions ri ,0 o 

It follows from the description of £l e and £lh that the degenerate region 
Oo is a union of strips parallel to the t axes and moreover every point of the 
boundary of each strip belongs, in fact to f^oo- Then, it follows that Oo — ^oo 
is an open set. We claim next that the degeneration is maximal everywhere, 
i.e. f^o = Ooo- Indeed, consider an integral curve of the A3-family lying 
inside this open set f^o — fioo- Then r% has constant values along the curve, 
and therefore it can not approach the boundary of the region VLq — fioo, 
since otherwise inside one would get r\ = r 2 = r% inside the region, which 
contradicts the assumptions. Thus the whole characteristic stays inside the 
region f^o — Ooo- But then, exactly as above, in the proof of the Theorem 
13.11 we have that the derivative {r%) q must explode in a finite time unless it 
vanishes. And therefore, in the whole region f^o — Oooj r 3 is a constant and 
therefore £Iq = Oqo, and we are done. 
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6. Proof of the Lemma 12.11 

It is rather simple to prove (jaj) and (jbj) of the lemma. By the definitions 
of Aj and r% we have 

F p (\i) = 0, F(\i) = n 

Differentiate these two formulas with respect to r,-. Notice that the roots Aj 
depend on rj and also the coefficients U{ of the polynomial F. We get 

(11) F P p(Ai)(Ai) r . + F m (Xi) = 0, F r , (Ai) = 6 tj 
Then write the following identity 

(12) F ri (p)(p-X l ) = (u 1 ) r ,F p (p) 

which becomes clear, if one notice that on both sides there are polynomials of 
the same degree n with the same leading coefficient (u\) n , and both having 
Ai, . . . , A ra as the roots (the right hand side-just by definition, and the left 
hand side-by (fTTJO . 

Differentiate this identity with respect to p to obtain 



(p - Xi)F m + F ri = (ni) n Fj 



!>!> 



Substitute in this formula p = Aj one arrives to (jcj) of the lemma immediately. 
Substitute p = Xj and take into account (jlip . then one proves fb]). 

In order to prove (jaj) let me use the relation F r .(X{) = Sij of ([TTj) in order 
to conclude, that F ri is, in fact, identical to the i-th Lagrange interpolation 
polynomial: 

Differentiate this identity with respect to p at the point A^: 
F P n(\) = (k)p(Xi) = 2^ n~TY) 

Using (fTl]) we obtain 

F pri (Xi) __ 1 y^ 1 



(A 

r.... r ...A A ! 



Fpp Fp P {Xi) ^. (Xi — X s ) 



This gives the proof of (jaj) 

In order to derive (jd|) write (jb|) in the form 

C\?')rj = u ri/(^i — Aj) 

Following Gibbons, Tsarev ([7]) differentiate this formula with respect to r^ 
to get 



(13) (Xj)nr k 



U ri U rk 



(Xi — Xj) (Xi — Xj)(Xk — Xi)(Xk — Xj) 
Now change the order of the indices i and k in (|13|) to have 

/i<\ /\\ u r k rj u r k u rj 

(Xk — Xj) (Afc — Xj)(Xi — Xk)(Xi — Xj) 
Subtract now ([13]) and ([Hj) . One gets ([d]) of the lemma. 
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7. Concluding remarks and questions 

1. It would be very interesting to generalize the analysis presented here 
for 3 by 3 system to the case of higher orders. In general, there are much 
more cases of degenerations of eigenvalues for higher n. This makes the 
analysis much harder. However, the case of strictly Hyperbolic solutions 
seems to be tractable. 

2. Besides the question on periodic solutions for the quasi-linear sys- 
tem, there is interesting question to understand the global existence of smooth 
(not necessarily periodic in time) solutions, having Hyperbolic initial data. 
In other words assume, that we are given for t = the initial condition, say 
periodic functions in q Uo(q) = (uio(q), . . . ,u n o(q)), such that the eigenval- 
ues of the matrix A{Uq) are real and distinct. Now switch on the dynamics. 
The question is if there are solutions existing for all times, or blow up of the 
solutions can be established. 

3. One of the key ingredients of the proof in this paper was the observation 
that the (semi-)Hamiltonian property of the quasi-linear system allows one 
to push forward the Lax analysis of the blow up of solutions. It would be 
interesting to know, what can be said in this perspective for other reductions 
of the Benney chain. 

4. It was proved by Tsarev f[14j). that in principle, the solutions for 
the system can be obtained by the so called generalized hodograph method. 
This method however, is very implicit and therefore, it is not clear to me 
how it can be used in order to answer the question of long time existence. 
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